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We discuss the universal transport signatures near a zero-temperature continuous Mott transition 
between a Fermi liquid (FL) and a quantum spin liquid in two spatial dimensions. The correlation- 
driven transition occurs at fixed filling and involves fractionalization of the electron: upon entering 
the spin liquid, a Fermi surface of neutral spinons coupled to an internal gauge field emerges. We 
present a controlled calculation of the value of the zero temperature universal resistivity jump 
predicted to occur at the transition. More generally, the behavior of the universal scaling function 
that collapses the temperature and pressure dependent resistivity is derived, and is shown to bear 
a strong imprint of the emergent gauge fluctuations. We further predict a universal jump of the 
thermal conductivity across the Mott transition, which derives from the breaking of conformal 
invariance by the damped gauge field, and leads to a violation of the Wiedemann-Franz law in the 
quantum critical region. A connection to organic salts is made, where such a transition might occur. 
Finally, we present some transport results for the pure rotor O(N) CFT. 
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I. INTRODUCTION 

Despite decades of study, the Mott metal-insulator transition remains a central problem in quantum 
many body physics.^ In recent years attention has refocused on an old question: can the Mott tran- 
sition at T = be continuous? Usually the Mott insulating state also has magnetic long range order 
or in some cases broken lattice symmetry which doubles the unit cell. In such situations a continuous 
Mott transition between a symmetry unbroken metal and the Mott insulator requires not just the 
continuous onset of the broken symmetry but also the continuous destruction of the metallic Fermi 
surface. It is currently not clear theoretically if such a continuous quantum phase transition can ever 
occur. 

Considerable theoretical progresses! has been possible in situations in which the Mott insulator 
does not break any symmetries but rather is in a quantum spin liquid (SL) stateP The evolution 
from the Fermi liquid (FL) metal to the quantum spin liquid state is of interest because it provides 
an opportunity to understand the fundamental phenomenon of the Metal-Insulator Transition (MIT) 
without the complications of the interplay with the onset of broken symmetry. Such a transition 
has acquired experimental relevance with the discovery of quantum spin liquid Mott insulators near 
the Mott transition in a few different materials, notably the quasi-two-dimensional triangular lattice 
organic saltsPHUl re -(BEDT-TTF) 2 Cu2(CN) 3 and EtMe 3 Sb[Pd(dmit) 2 ]2, and the three dimensional 
hyper-Kagome material Na4lr 3 0s- Indeed upon application of hydrostatic pressure, these quantum 
spin liquids become metallicP3(The phase diagram of k-(BEDT-TTF)2Cu2(CN) 3 in addition has 
superconductivity at lower temperatures.) The nature of the transition in the experiments is not 
currently understood but could be potentially described by a SL-FL quantum critical MIT. 

A wide variety of quantum spin liquid phases can exist theoretically. However a natural stat d 3 * 18 ^that 
emerges near the Mott transition is a gap less quantum spin liquid which has a Fermi surface of spin- 1/2 
neutral quasiparticles, the spinons, while the charge excitations are fully gapped. Consequently there 
is a single particle gap in the electron spectral function even though there are gapless spin-carrying 
excitations. This type of spin-charge separation where only the charge localizes can be favored in 
systems (such as the organics) where frustration and charge fluctuations suppress magnetic ordering. 

The MIT we consider was studied at the mean field level in Ref. A subsequent analysis^ of 
the quantum fluctuations provided evidence that the second order nature can survive the (inevitable) 
inclusion of many-body effects. A rich set of properties associated with the quantum critical point 
(QCP) was uncovered, many not present at the mean field level. For example, it was found that 
quantum fluctuations lead to a divergence of the effective FL-mass as one approaches the MIT. This 
gives rise to a divergence of the specific heat capacity coefficient, 7 = C/T. A key insight was the 
realization that the quantum critical state at the edge between the FL and SL is actually a non-FL 
(11FL) metal where the Landau quasiparticle is destr oye d but the concept of a sharp Fermi surface 
persists. This was dubbed a "critical Fermi surface" PEI One prediction of this theory was that the 
zero-temperature resistivity jumps by a finite amount as one goes from the FL to the quantum critical 
state, where the value of the jump was predicted to be universal, Rh/e 2 , R being a dimensionless 
number associated with the QCP. This resistivity jump is illustrated in Fig. [TJ 

The main purpose of this work is to provide a controlled calculation of the value of this jump 
and, more generally, analyze the behavior of the electric resistivity in response to changes in pressure 
(modifying the metallic bandwidth) and/or temperature in the vicinity of the QCP. The resulting 
resistivity can be compared with experiments, where the bandwidth can be changed by applying me- 
chanical or chemical pressure. The predictions we make, viz. sharp resistivity variation on the order of 
I0h/e 2 r~j 100 kO, quantum critical collapse of pressure and temperature dependence, thermal conduc- 
tivity jump and violation of Wiedemann-Franz law, provide distinct signatures of the fractionalization 
at the MIT and of the critical Fermi surface state intervening between the FL and SL. 

A crucial ingredient of the theory is that, at zero temperature, the emergent gauge fluctuations 
associated with the fractionalization decouple from the quantum critical charge fluctuations. Despite 
this the gauge fluctuations are expected to play a crucial role for non-zero temperature transport 
properties. A similar phenomenon happens in the Kondo breakdown model studied in Ref. 1191 An 
important difference with the Kondo breakdown scenario is that the charge fluctuations near the 
Mott transition studied in this paper are described by an interacting field theory at low energies. 
Irrespective of this difference, the same conclusion holds: the gauge fluctuations become important 
for low frequency transport at non-zero temperature. We show this explicitly by calculating the effects 
of these gauge fluctuations on the transport. In particular the precise value of the universal resistivity 
jump in the limit that the frequency of the applied electric field goes to zero faster than temperature 
is strongly affected by the gauge fluctuations. In contrast, in the opposite order of limits the universal 
resistivity jump is unaffected by the gauge field. We further predict a universal jump of the thermal 
conductivity across the Mott transition, which derives from the breaking of low energy conformal 
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FIG. 1. Jump of the universal resistivity, p, and Lorentz number, k/ctT, at T = as a function of 8, 
which is proportional to ratio of the bandwidth to the Hubbard repulsion, Eq. The latter jump signals a 
violation of the Wiedemann- Franz law by the critical Fermi surface state, k is the thermal conductivity; R, K 
are universal constants associated with the Mott QCP. In particular, they strongly depend on the emergent 
gauge boson associated with the electron fractionalization. The resistivity becomes infinite in the SL, and as 
a consequence so does the Lorentz number. 
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FIG. 2. The phase diagram of the quantum critical Mott transition. S tunes the ratio of onsite repulsion 
to the bandwidth away from its critical value, and can be put in correspondence with P — P c , the deviation 
from the QC pressure, P c . The dark shaded (blue) region is the quantum critical region, where the Landau 
quasiparticle is destroyed but a "critical Fermi surface" nonetheless exists. It separates the spin liquid (SL) 
and the Fermi liquid (FL). The intermediate-T states (with prefix 'M'), the marginal SL and FL, differ from 
the low temperature ones by the fact that the spinons and gauge bosons still behave as in the QC region. 



invariance by the gauge field, and leads to a violation of the Wiedemann-Franz law by the critical 
Fermi surface. 

The paper is organized as follows. In section [TlJ we summarize our main findings; section |III| 



introduces the slave-rotor description for the Hubbard model. In section IV we formulate the transport 
via a quantum Boltzmann equation for the critical charge fluctuations, the rotors, and present its 
solution at criticality. Section [V] extends the resistivity calculation to the entire QC region. In |VI| 
we discuss the behaviour of the resistivity at large temperatures and frequencies. Signatures relating 
to thermal transport are discussed in section |VII[ The Appendices give further details regarding the 
critical rotors, with a focus on their transport properties in the absence of the emergent gauge field, 
i.e. in the pure O(N) non-linear sigma-model (NLctM) for N > 2. 



II. MAIN RESULTS 



In this section we summarize our main results. We first need to briefly describe the finite temper- 
ature phase diagram obtained in Ref. [5J which is reproduced in Fig. [2j The parameter S tunes the 
ratio of the electronic bandwidth to the onsite repulsion away from its critical value, and can be put 
in correspondence with P — P c , the deviation from the quantum critical pressure, P c : 



8<xt/U-(t/U) c 



P-Pr 



(1) 



At T > 0, the metal-insulator transition becomes a crossover due to the presence of the emergent 
gauge boson. We distinguish three main phases in Fig. [2j the spin liquid (in white in the figure), 



4 



the Fermi liquid (pink/light shading), and the quantum critical state bridging the two (blue/dark 
shading). The latter state is a non-FL where the Landau quasiparticle has been destroyed, yet a 
sharp Fermi surface persists: an instance of a "critical Fermi surface" . In exiting the QC region, one 
enters two intermediate phases: a marginal spinon liquid (MSL) or a marginal Fermi liquid (MFL). 
These are similar to their low temperature counterparts, the SL and FL, except that the spinons and 
gauge bosons still behave as in the QC region. As these correspond to fluctuations in the spin degrees 
of freedom, the two crossovers may be interpreted as corresponding to spin and charge degrees of 
freedom exiting criticality at parametrically different temperatures. At sufficiently low temperature 
they crossover to the usual SL and FL states. The behavior of the electric resistivity as one tunes 
across the phase diagram is illustrated in Fig. [3] Panels a) and c) correspond to the T-dependent 
behaviour at fixed S (i.e. pressure), and vice- versa for b) and d). The important crossovers for low 
temperature transport are the boundaries of the QC region: there, the charge degree of freedom either 
localizes (SL) or condenses (FL). At the former crossover the resistivity becomes thermally a ctiva ted, 
~ e A +/ T , because of the finite Mott charge gap A + . This can be seen in curves 1-2 in Fig. 3(c) At 
the crossover to the FL, it abruptly drops to its residual metallic value p m (curves 4-5 in Fig. 3(c) I. 
The regime of interest for transport corresponds to the QC non-FL, where the resistivity relative to 
its residual value in the metal, p m , is purely universal: p — p m « (h/e 2 )R, where R is a universal 
dimensionless constant. Our controlled calculation of R in a large- A approximation gives the estimate 
R = 49.8. R sets the size of the jump shown in Fig. [lj which is reproduced in Fig. |3(d)l curv e 1 



At finite temperature, this jump becomes a steep increase, as shown in curves 2-3 of Fig. |3(d) We 



emphasize that the low temperature resistivity above the QCP, S = 0, is T-independent and takes the 
value p = p m + (h/e 2 )R. 

The diverse behavior shown in Fig. [3] can be obtained from a single- variable function. Indeed, the 
temperature and pressure dependent resistivity (relative to its constant residual value in the FL) can 
be collapsed by a universal scaling function G associated with the Mott QCP: 



= iG ( — ) , (2) 



where the dynamical and correlation length exponents correspond to those of the 3D XY universality 
class: z = 1 and v s» 0.672. Indeed, the critical charge degrees of freedom can be effectively described 
by a Bose-Hubbard model at half-filling near its insulator-superfluid transition, which belongs to that 
universality class. We show that although the spin fluctuations encoded in the emergent gauge field 
associated with the electron fractionalization do not alter these exponents, they have strong effects 
on the scaling function, and thus on the value of the universal jump, (h/e 2 )R. 

We predict that thermal transport also shows signatures of the critical Fermi surface. In particular, 
the thermal conductivity divided by temperature, n/T, has a universal jump at criticalit y, by an 
amount (kg/h)K, where K is a. dimensionless number just like R. As we explain in section VII| the 



emergent gauge fluctuations play an important role by breaking the conformal invariance present in 
their absence, thus reducing k/T from a formally infinite value to a finite, universal one. Finally, 
combining the electric resistivity and thermal conductivity jumps, we predict that the QC non-FL 
violates the Wiedemann-Franz law by a universal amount: the Lorentz number differs from its usual 
value in the FL by (/cs/e) 2 i?AT, as shown in Fig. [I] 



III. MOTT TRANSITION IN THE HUBBARD MODEL: A SLAVE-ROTOR 

FORMULATION 



To set the stage we briefly review the description of the insulating quantum spin liquid with a spinon 
Fermi surface,^ and the continuous bandwidth-tuned Mott transition^ to it from a Fermi liquid. We 
consider a single-band Hubbard model at half-filling on a 2D non-bipartite lattice (for e.g. triangular): 

H = -tJ2 (4r<W' + h.c.) + Uj2( n r - I) 2 , (3) 

(rr') r 

where c ar annihilates an electron with spin a at site r, and n r = cl r c ar . In the small U/t limit, 
the ground state is a Fermi liquid metal, while in the opposite limit a Mott insulator results. The 
interplay of frustration and strong charge fluctuations can lead to a quantum spin liquid ground state 
instead of a conventional antiferromagnetic Mott insulator. We shall focus on the transition to such 
a state. 

The slave-rotor construction^ is tailor-made to describe the spin-charge separation that occurs as 
the charge localizes when the electronic repulsion becomes sufficiently large, yet weak enough for the 
spins to remain disordered, even at T — 0. At the level of the microscopic Hubbard model, Eq. pi), 
the slave-rotor construction is a change of variables to degrees of freedom better suited to describe 
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FIG. 3. Sketch of low temperature behaviour of the resistivity near the quantum critical (QC) Mott transition. 
Panel c) shows the resistivity vs T for different values of the onsite repulsion over the bandwidth (tuned by 5), 
with the corresponding cuts shown in the phase diagram in a). Panel d) shows the resistivity vs S at different 
temperatures, with the corresponding cuts shown in the phase diagram in b). In c) and d), the markers 
correspond to the location of the resistivity jump upon entering the QC state from the FL. The value of the 
jump is universal: Rh/e 2 . Our calculations yield R = 49.8, which translates to a jump of ~ Sh/e 2 . 



the SL, in which the electron is fractionalized into spin and charge carrying "partons" : 

Crcr = IpvrW ■ (4) 

The fermionic spinons, i/w, carry the spin, while the bosonic rotors, b r — e~ l6r , the charge of the 
original electron. The projection from the enlarged Hilbert space to the physical one is obtained from 
the operator identity relating the rotor charge or "angular momentum", to the fermion number, nf. 
lb = 1 — nf, which is enforced at each site, where nf = n is the actual electronic occupation number 
(because \b\ = 1). By virtue of Pauli exclusion, the charge relative to half-filling at each site can 
only be —1 (double occupancy), +1 (hole) and (single occupancy). Hence, the positive (holon) and 
negative (doublon) charge excitations encoded in the rotors relate to the holes and doubly-occupied 
sites of the half-filled Hubbard model, see Fig. |4j Moreover, since the system is at half-filling, there is 
a low-energy particle-hole symmetry between these positive and negative charge excitations. 

In the long-wavelength limit, a U(l) gauge structure emergesP The temporal component of the 
gauge field results from the above constraint necessary to recover the physical Hilbert space, while 
the spatial components derive from the fluctuations of spinon bilinears about their saddle-point con- 
figuration. After coarse-graining, the low-energy effective action for the Hubbard model in terms of 
the fractionalized degrees of freedom can be written as 

S = Sf,M + Sf^ a + S a , (5) 
S b . a = (K d » ~ m -) h \ 2 + iA d & | 2 ~ !)) > ( 6 ) 

Sf.a = J Vv (d T - /i - ia + ^2m <1 " > ) ^ a ' ^ 
S a = \ I (e^d 7 a p ) 2 . (8) 

e Jx 

We work in units where the rotor velocity, c, is set to one, unless otherwise specified. The complex 
boson field b is constrained to lie on the unit circle via the Lagrange multiplier field A. The indices 
v, 7, /3 run over imaginary time and the two spatial dimensions; \i is the electronic chemical potential. 
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We have used the shorthand f — L dr J d 2 x. The parameter that tunes the Mott transition of the 
rotors (and hence of the whole electronic liquid) is g oc U/t, where U is the Hubbard repulsion from 
the original electronic Hamiltonian, while t is proportional to the electronic bandwidth. To make the 
action dimensionless, the parameter g carries the dimension of length, as given by a real space cutoff 
scale. We can relate it to the parameter 5 introduced in Eq. ([!]) via 5 = g^ 1 — g^ 1 oc t/U — (t/U) c . 
For small coupling, g < g c , the rotors spontaneously condense, corresponding to the metallic phase 
of the original Hamiltonian; while in the opposite limit, g > g c , the rotor field is disordered leaving 
the system in a SL ground state. The gauge fluctuations have a Maxwellian action resulting from the 
elimination of high energy fluctuations; eg is the corresponding bare gauge charge. 

The emergence of a relativistic action for the rotors, which have a dynamical exponent z = 1, is a 
consequence of the emergent low-energy particle-hole symmetry of the Hubbard model at half-filling 
noted above. This low-energy symmetry will be important when we examine the effect of the gauge 
field on the critical rotors. It will lead to a strong suppression of the dynamical gauge fluctuations in 
the charge (rotor) sector. 

The field theory above is strongly interacting. Indeed, in two dimensions the rotor NLuM considered 
separately flows to a strong coupling fixed point where even the 6-field quasiparticles are ill-defined. 
The spinons and gauge fluctuations do not alter this. One perturbative approach to the problem 
extends the field theory to include a large number of flavors of the matter fields: when that number 
is very large, we have weakly interacting quasiparticles, at least in the boson sector. We will use this 
extension, which we now describe in more detail, to bring the calculation under control. 



A. Low energy theory and large-iV extension 



We consider the slave-rotor field theory extended to have a large number of rotor and spinon flavors, 
allowing for a systematic study of transport The number of copies of the complex rotor is taken to 
be N/2, yielding an even number, N, of real scalar fields. The case of physical interest has N = 2. In 
this large- AT extension, the effective actions of the rotors and spinons read: 

N f 

S b , a = Y g j (\( d » - ia ») b «\ 2 + ZA (M 2 - !)) . ( 9 ) 

Sf, a = Jj>„ [dr -li- ia a + (V 2to ? " )2 ) Vv , (10) 

where each b a is a complex scalar, and a runs from 1 to N/2, while there are N copies of spinons, 
a = 1, ...,N. Repeated indices are summed over, for example \b a \ 2 = b a b a = J2 a \b a \ 2 - A is the 
Lagrange multiplier field enforcing the constraint that the O(N) real field be unimodular. The coupling 
g has been rescaled g —> g/N. Note that the gauge field reduces the rotor global symmetry from 0(N) 
to U(N/2). In the N — > oo limit, the fluctuations of the A and gauge bosons are unimportant. Formally 
integrating out the rotors yields the following form for the partition function 20 Z — J V\e~ s ° !! ^ , with 



trln(-<9 2 

9 



(11) 



The overall factor of N plays the role of 1/h so that in the large- N limit the quantum fluctuations 
are suppressed and we only need to consider the classical equation of motiorPSl 



9 1 - =-, (12) 
q z + m z g 



where J — J d d q/(2ir) d and the mass squared, m 2 = iAo G R, corresponds to the saddle point 

value of the uniform component of A, Ao- This mean- field value of A plays the role of the mass for 
the rotors in their insulating phase at large g. It can be alternatively seen as the inverse correlation 



length, £ ~ 1/m. At sufficiently small g, Eq. (12) has no solution, and a different approach must be 



used to describe the condensation. In the following, we shall focus mainly on the quantum critical 



regime as well as on the insulating phase. The solution of the saddle-point equation, Eq. (12), in 



the large- N limit and at finite temperature directly above the quantum critical point, g — g c , yields 

m = GT, where 8 = 2 ln( ~ 0.96, twice the logarithm of the golden ratiopHThe mass vanishes 
linearly as T — > 0, i.e. the correlation length of the charged rotors diverges upon approaching the 
QCP: £ ~ 1/T. The full dependence of m on g and T at N = oo is given in section |vj 

Corrections at order 1/N to the N — oo saddle-point correspond to interactions mediated by the 
A and a bosons, which develop dynamics when N is finite. The rotor action, including the effective 
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mass corresponding to the saddle-point value of the A boson, now reads 



i v )K\ 



m 2 \b a \ 2 +i\\b a \ 



(13) 



We are using the Coulomb or transverse gauge so that it is understood that we only include con- 
figurations where V ■ a = 0. The transverse and temporal component of are decoupled in this 
gauge. Further, we can omit the latter in the low-energy limit because it is screened by the spinon 
Fermi surface. The transverse part of the gauge field remains gapless because the currents remain 
unscreened, as opposed to the charge. In the remainder, we shall use a to represent the transverse 
component. The 1/iV corrections to the saddle point will generate 0(1/ N) propagators for both the 
gauge and A bosons, which acquire the following effective action: 



Hq)\ 



,N 



\Hq)\ 



-IT 



2 \ ' "J 1 ^" 2 
where the finite temperature, imaginary-time polarization functions read 

d 2 p 1 1 



(2tt) 2 K 
d 2 p 



vi) 



(2q X p) 5 



1 



(2tt) 2 K + ^) 2 +£p+ g f% + 4 



(14) 

(15) 
(16) 



where the superscript "j" identifies the current-current correlator; we have defined the rotor dispersion 
relation: e p = ^/p 2 + to 2 . Details about the computation of lib, n J b can be found in Appendices |a| and 
[c| respectively. II J b is discussed in the next subsection, IIIAl 



n 



fc 2 

k 771 



The spinon Fermi surface contributes 
,| < v F k (17) 



where the Ci are real numbers, while kp,VF are the Fermi momentum and velocity, respectively. As 
we work in units where the velocity of the rotors is set to 1, we need to keep vp explicitly. Because 
of the term \ui n \/vFk in the fermionic polarizability, the gauge fluctuations are Landau damped. 



1. Role of gauge fluctuations 

We now examine the role of the gauge fluctuations. The Landau damped dynamics due to the 
Fermi surface dominate those induced by the rotor polarization function, n J fc , and we can evaluate the 
latter in the static limitP Note that Hi, just as IT, depends on the temperature and g via the mass 
of the rotors, to. For g = g c , as shown in Appendix [Cj we get 

[a%q if<7>T 

where 72 ~ 0.031 is a dimensionless constant, while cr" w 0.063 is the rotor condu cti vity in units of 
e 2 /h in the large frequency or T — > limit w/T» 1. As discussed in sections IV A and VIA it differs 




from the DC conductivity we are seeking and can be obtained from a simple T = analysis. The 
q 2 behaviour at q -C T results from including the mass of the rotors, to = 0T, when computing the 
current polarization function. It has the behaviour expected from massive modes since for q <C T, 
the fluctuations exceed the correlation length, 1/m ~ 1/r, and must be gapped. The important 
term for the low-temperature behaviour is the linear q contribution. This non-analytic dependence 
arises because of the gaplessness of the critical rotors and gives the gauge field a, z a = 2 dynamical 
exponent, making it less singular than deep in the SL where the rotors are gapped and the gauge 
bosons have z a — 3. Indeed, the z a — 2 damped gauge fluctuations give the spinons a self-energy 
~ jjiu)ln(/j,/\u}\), which is weaker than in the usual SL, where we have -^i|u;| 2 / 3 . Thus, at criticality, 
as well as in the MFL and MSL phases, the spinons form a marginal Fermi liquid, leading to the usual 
logarithmic corrections. We emphasize that we do not need to worry about possible subtleties with the 
breakdown of the naive large- N expansion for a Fermi surface coupled to a gapless bosom 2 ^ Indeed, 
we are mainly concerned with the quantum critical region, where, as stated above, the fermions only 
acquire logarithmic corrections due to gauge fluctuations. Such a marginal Fermi liquid of spinons can 
be controlled by a simple perturbative renormalization group (RG) approach! 23 ! 24 ! In contrast, deep 
in the SL, one might need to take the limit of small z a — 2 simultaneously with 1/N — > to make 
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(a) (b) 



FIG. 4. Charge excitations near the Mott transition, a) Triangular lattice at half-filling; the small shaded 
disks represent electrons. The double-occupied (empty) sites are identified by a red/left (blue/right) circle. 
These are encoded in the charge rotor excitations, the holons and doublons, respectively. Under an applied 
electric field, they will move in opposite directions, b) By virtue of the emergent particle-hole symmetry 
between doublons and holons, it is possible to have a state with zero momentum, P, but finite current, J. 
This allows interactions to dissipate current while conserving momentum. 

the expansion controlled.^ Further the main transport properties will derive from the rotor sector for 
which the large- N works reliably. The spinons affect the rotors only via the damping of the gauge 
bosons, which we believe is a robust feature, independent of the expansion scheme. 

Regarding the rotor or charge sector, it was showrP that the gauge fluctuations do not alter the 
nature of the rotor excitations, i.e. the rotor self-energy is subleading compared to the bare dynamics. 
Insofar as the thermodynamic critical properties of the charge sector are concerned, they belong to the 
3D XY universality class, unaffected by the gauge bosons or spinons. The importance of the damping 
at quenching the gauge fluctuations can be heuristically understood by examining the dominant rotor 
fluctuations, which have u> *~ q (z = 1). Substituting this dispersion relation into the Landau damping 
term we get p\ui\/q ~ p. Thus, the dominant rotor fluctuations see the gauge bosons as screened. 
Such an effect was also identified in Ref. 1251 in the context of a quantum critical transition between a 
Neel-ordered Fermi-pocket metal to a non-FL algebraic charge liquid, called a "doublon metal" . This 
suppression mechanism of the gauge field due to a Fermi surface was referred to as a "fcrmionic Higgs 
effect". 

In the next section, we shall show that although the gauge fluctuations are not effective at influencing 
the thermodynamic critical properties in the charge sector, they can have strong effects on non-zero 
temperature transport, down to arbitrarily low temperatures. 



IV. CRITICAL TRANSPORT NEAR THE MOTT TRANSITION 

In a slave-particle theory such as the one under consideration, many observables can be determined 
from the separate responses of the partons. These relations generally go under the name of Ioffe-Larkin 
composition rules. For example, the one for the resistivity reada^ 

P = Pb + Pf, (19) 

where pb./ is the resistivity of the spinons and rotors, respectively. The resistivities "add in series" 
because of the constraint relating the spinons and rotors to recover the original Hamiltonian/Hilbert 
space: The electric field induces a motion of the electrically charged rotors, forcing spinons to flow as 
well. Alternatively, we can say that the external electric field induces an internal one. It follows that 
the parton with the highest resistivity governs the entire electric response. Near the Mott transition, 
the rotors have the most singular response as they undergo a quantum phase transition, while the 
spinons form a Fermi surface throughout. We thus anticipate that the strong variation of pb across 
the transition will give the entire resistivity its key dependence. 

Let us first discuss the T = and clean limit, in which case the spinon Fermi surface has vanishing 
resistivity throughout. The rotors also have vanishing resistivity in their condensed phase, such that 
the FL has j o = pt + /?/=0 + = 0as expected. On the Mott side, the rotors are gapped hence 
the whole system has infinite resistivity: p\, — oo = p. The interesting feature happens directly at 
criticality, where altho ugh w e still have pf = 0, the rotors have a finite universal resistivity induced 
purely by interactions j 21 l 27 l p b = Rh/e 2 . It is possible for systems with particle-hole symmetry or 
cquivalcntly emergent relativistic invariance to have a finite resistivity in the absence of disorder 
or Umklapp scattering. For such systems, the momentum and electric current operators need not be 
proportional, allowing interactions to dissipate the latter while preserving the former. More physically, 
these systems have independent and symmetry-related positive and negative charge excitations that 
flow in opposite directions under an applied electric field yielding a state with a finite current but 
with zero momentum. This is schematically illustrated in Fig. [4] The finite, interaction-driven rotor 
resistivity at criticality leads to a discrete jump at T = 0, as illustrated in Fig. [T] 
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This scenario naturally extends to finite but low temperatures: the fermions still contribute only 
a constant pf, which is zero for a clean system, or finite in the presence of weak disorder. Instead 
of discontinuously jumping, the rotor resistivity increases rapidly upon entering the quantum critical 
region, where it slowly increases until the growth becomes exponential at the crossover to the spin 
liquid, as is shown in Fig. |3(d)| We shall thus focus on the resistivity of the rotors, pb, for which we 
perform a 1/N expansion. In the simplest limit, TV = oo, the rotors are free because they decouple 
from the A and gauge bosons. The DC resistivity thus vanishes in the absence of scattering, pb = 0. 
At order 1/N, the rotors begin colliding with the constraint field A and the emergent gauge boson, 
leading to a finite resistivity. For sufficiently large TV, the system has well-defined quasiparticles, whose 
transport properties can be unambiguously studied by a quantum kinetic (or Boltzmann) equation, 
to which we now turn. 



A. Quantum Boltzmann equation for critical charge fluctuations 



We formulate the quantum Boltzmann equation (QBE) for the distribution functions of the rotor ex- 
citations in the presence of an oscillating electric field, E(t), with driving frequency uj. This frequency 
plays an important role as it introduces an energy scale that divides the critical frequency-dependent 
resistivity into two regimes: ui < T and w > T. Indeed, the frequency w must be compared with 
the dominant scale for the rotors in the quantum critical regime, which is the temperature. As was 
established in seminal work by Damle and SachdevpD the limits u — > and T — > do not generally 
commute for the response functions of critical systems. For example, the T = DC resistivity is ob- 
tained by first taking uj/T — > 0, then T — > 0, so that one must necessarily perform a finite temperature 
analysis to obtain the correct DC response. A T = calculation, which is equivalent to taking T — » 
first, yields the response in the u)/T — > oo limit, which generically differs from the DC behaviour. This 
non-commutativity of the uj — > and T — > limits can be explained on physical grounds: The small 
frequency resistivity (uj < T) is dominated by the incoherent scattering of thermally excited critical 
fluctuations; it corresponds to the hydrodynamic limit. In contrast, the large frequency resistivity 
(uj > T) arises from the coherent motion of field generated excitations; it is mainly collisionless. The 
dichotomy is even more striking in our case due to the presence of the gauge bosons: We shall show 
that although the gauge fluctuations do not affect the transport in the large uj/T limit, they actually 
dominate the DC resistivity! 

We assign the electric charge to a single rotor flavor, b±, which couples to the oscillating electric 
field E(t). The standard mode expansion for the electrically-charged rotor operator reads: 

&x (a?) = I a + (t, k)e ik x + al (t, k) e - %kx , (20) 

Jk 

where we have defined a±/a± as the annihilation/creation operators for holons (+) and doublons 
(— ), i.e. the positive and negative electric charge excitations. The expectation value of the current 
can be decomposed into two pieces: J(t) = Ji(t) + Jn(t), where 

J 'W- / y>-<at(t,fc)a a (t,fc)) (21) 
= /]>>-/ s (i,fc) (22) 

J k „ l e k 



s=± 



Jn(t) = [ ^-{aUt,-k)aUt,k) - al(t,-k)al{t,k) - 2aUt,-k)al(t,k)) + h.c. (23) 
Jk 2e fe 

We have defined the distribution functions of positive and negative charge e xcit ations: /. 



{cv s (t, k)a s (t, k)), s = ±; efc = V m 2 + k 2 is the rotor dispersion. From Eq. (23), it should be 
apparent that as Jn involves pair production, it will only contribute when the driving frequency is 
above the pair-production threshold, uj > 2m, where m ~ T in the QC region. We shall concern our- 
selves with the determination of «//, i.e. f s (t, k), which governs the transport in the small frequency 
limit. The asymptotic high frequency resistivity in the limit uj 3> T can be obtained from a T = 
calculation and we leave its analysis to section |VI A| 



The QBE for the distribution function of holon (doublon) rotor excitations, f s with s = ±, respec- 
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FIG. 5. Illustration of the main scattering mechanisms determining the resistivity in the QC region. The blue 
disk corresponds to a holon excitation with charge +e. In addition to the usual scattering between critical 
charge fluctuations (mediated by the A field), the static emergent gauge fluctuations generate a random 
"magnetic field", V x a, that scatters the holons and doublons. A schematic configuration of this emergent 
magnetic field (which is always perpendicular to the plane) is shown, where the scale gives its strength and 
direction, the latter dictated by the sign. 



tively, reads 



(d t + SE ■ d k )f s (k,t) = -(h[f±] + Ia[f±] 



dn f d 2 q 



it J (2tt) : 
(2ir)6(e k - e k+q - O) 



u b (n, q ) 



(2k x qYQD(n,q) 



(24) 
(25) 



+ 



4e fe e fe+9 
(2n)S(e k - e k+q + Q) 

(2n)8(e k + e_ k+q - Q) 



f s (k,t)(l + f s (k + q,t))(l + n(Q)) 

f s (k + q,t)(l + f s (k,t))n(n) 
f s (k,t)(l + f s (k + q,t))n(n) 
f s (k + q,t)(l + f s (k,t))(l + n(Q)) 
f s (k,t)f. s (-k + q,t)(l + n(n)) 



(1 + f s (k, i))(l + f_ s (-k + q,t))n(Sl) 



We have absorbed the magnitude of the rotor charge into E. The right hand side of the QBE, the 
collision term, can be obtained for instance by invoking Fermi's Golden RulelMSH xhe £ rs ^ ^ wo §. 
functions enforce energy conservation for absorption and emission of A and gauge bosons by the rotors, 
while the last one corresponds to pair creation/annihilation of holon-doublon pairs. 

The propagators of the A and gauge bosons, 2/7VIT b and 2D/N = 2/iV(ffi +n J b ), respectively, enter 
into the QBE via their spectral functions which dictate the density of states the rotor excitations 
can scatter into. They are evaluated in equilibrium. This is justified in the large-A^ limit since the 
external field couples to a single rotor flavor such that the associated non-equilibrium corrections to 
the polarization functions sublead in \/N . In other words, the drag of the A and gauge fields by the 
electric field, an analogue of "phonon drag", is negligible in the large- N regime. The rotor flavors 
that do not directly couple to the electric field, b a> i, play the role of an effective bath at equilibrium, 
from which the constraint field acquires its dynamics. 



The scattering terms on the RHS of Eq. ( 24 1 all scale like 1 /N because the gauge and A bosons 
have propagators of that order. As N — > oo, the scattering terms vanish and the rotors become free, 
displaying a sharp Drude peak in the real part of the small frequency conductivity: a b oc 8(uj), lj < T. 
As we shall see, the finite 1/N effects will cure this singularity yielding a finite DC conductivity. 

We now proceed to the solution of the rotor QBE, Eq. (24), by first expanding the distribution 
function to linear order in E, 



f s (k, bS) — n(e k )2TrS(uj) + sE ■ kip(k, oj) . 



(26) 



where we have Fourier transformed from time to frequency. The deviation function, ip(k,u>), only 
depends on the magnitude of k, since E ■ k fully encodes the rotational symmetry breaking in the 
presence of the external electric field. The unknown function (p parametrizes departures from the 
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equilibrium Bose-Einstein distribution, n(efc), with e 2 = m 2 + k 2 . The result of the linearization 
of the collision term due to the A bosons can be found in Appendix [X] we cannot make significant 
simplifications there. Indeed, we need to perform a careful numerical evaluation of the retarded 
polarization function ILj(f2, q) for all frequencies and momenta. On the other hand, we can significantly 
simplify the scattering term due to gauge bosons. As mentioned above, dynamical gauge fluctuations 
are suppressed in the rotor sector at T = 0. Although a static (CI = 0) gauge mode, a(0,q), escapes 
the Landau damping, at T = it constitutes a set of measure zero in the continuum of excitations 
and is thus unimportant. At finite T, the static Cl n = Matsubara frequency is well separated from 
the others and thus provides a viable scattering channel. As long as the hierarchy to <C T <C \i is 
maintained, where fi is the electronic chemical potential, this effect remains. We can interpret the 
situation as follows: the rotors are scattered by a static random magnetic field, V x a(x), generated 
by the emergent gauge fluctuations, which increases the resistivity compared with the usual insulator- 
superfluid transition of rotors. The width of this random distribution of static magnetic fields is 
proportional to the temperature T. This is illustrated in Fig. [5j To determine the corresponding 
scattering rate, let us first rewrite the gauge-boson-rotor scattering term: 



/ (27r)«5(e fe - e k+q - CI) 



(27) 



(27r)i$(e fc - e k+q + CI) 

&k e k+q 



f s {k,t){l + f s {k + q,t)){l + n{Cl)) 

f s (k + q,t)(l + f s (k,t))n(Cl) 
f s (k,t)(l + f s (k + q,t))n(Cl) 



-f s (k + q,t)(l + f s (k,t))(l+n(Cl)) 



We have neglected the particle-antiparticle production term, the one with 6(ek + e-fe+g — CY), because 
it requires an energy CI > 2m, which renders it subleading due to the suppression of dynamical gauge 
fluctuations. Linearizing yields the simple relaxational form 



sE ■ kip(k, uj) 



T a (k) 



Sfs 

"r»(A) 



with the momentum-dependent scattering rate 

1 T I 



,2 „,2 

dy 



r a (k) N ne k J n J b (0,2fcy) 



y V 1 - v 2 



(28) 



(29) 



This result for the scattering rate due to the emergent gauge bosons is valid for temperatures T <C fi. 
This elastic scattering rate obtained in the static regime is universal in the sense that it does not 
depend on the Fermi surface information, such as kp and vp. It follows that we can express it in 
terms of a single-parameter scaling function: 



1 



T ( k 
N Fa \T 



(30) 



The exa ct momentum dependence can be determined numerically and is shown in Fig. [6} In evaluating 

= result II J b (0,<7) ~ q, the rate would have vanished as 



Eq. (29|), had one simply used the T 



After linearizing Eq. (24), we get the following equation for the deviation, 99, from the equilibrium 
distribution: 



^(p,u)+g(p)/T 2 = ^ { - [F\M» + F„U»]^i>.^)+ / ./// l<xU>-i> > P U> ■ - ' 



(31) 



where we have Fourier transformed from time to frequency and have introduced dimensionless mo- 
mentum variables p,p' that have been normalized by temperature: p — k/T. On the LHS, the 
term that makes the equation non-homogeneous reads g(p) — d ep n(e p )/e p = — e ep /e p (e ep — l) 2 , 
where it is understood that when we use the dimensionless momentum, the rotor mass is scaled 
by T, e p = ^p 2 + (m/T) 2 , and the Bose function does not contain the usual factor of temperature: 
n(e p ) — l/(e £p — 1). We have also introduced the dimensionless kernel K\(p,p') describing the non- 
elastic processes in which the rotors exchange energy with the A bosons. It is independent of the 
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k/T 



FIG. 6. The scattering rate due to gauge fluctuations (solid line), Eq. (29 1, and due to the A-bosons, 
(dashed). The latter, Fx, appears in Eq. (311 and is defined in Eq. (A3 1. We see that the rate due to the 
gauge fluctuations is a factor of ~ 8 larger than that due to the A field and is thus the dominant scattering 
mechanism. 



gauge field, the driving frequency and N. In analogy with Eq. (30), we have further defined the 
dimensionless function Fx that corresponds to the scattering rate due to interactions with the A field. 
More details about F\ and K\ can be found in Appendix |A"} 



B. Solution of QBE and rotor conductivity 



By performing the rescalings 

*=f, (32) 
$(p,w) = ^(p,w), (33) 
where again p = k/T is the dimensionless momentum, we obtain a universal, parameter- free equation: 
-iQ$(p,L>)+g(p) = -F(p)$(p,Q) + J dp'K x (p,p')$(p',6j). (34) 

The gauge fluctuations do not spoil the existence of such a universal equation, which arises for the 
pure rotor theory J22l because they contribute a universal scattering rate modifying Fa — > F\ + F a =: F. 

The above integral equation needs to be solved numerically. Once we obtain $, we can compute 
the low- frequency conductivity, tif,/, from the expression for the current in terms of the distribution 
function of the rotors 



(35) 



= E*/(0^-M^) (36) 

Assuming the E field is in the x-direction, we get: 

a bI (oj) = (J Jx (u)) /E x (u) (37) 

l/VVC*^ (38) 
z7r Jq eh 

= iVxi/ A %?^ (39) 



2tt 

where A is the momentum cutoff used in the numerical solution. The last equality makes use of the 



13 



scaling function for ip, so that the small-frequency conductivity, ui/T <C 1, can be written as 

, U H = ^XE;S , (40) 

where the fundamental constants e and h were re-introduced; S/ is a complex valued universal function 
defined by Eq. ( 39 ) , the real part of which is shown in Fig. [7] As shown there, the conductivity 



is substantially reduced at small frequencies uj compared with the pure 0(N) model due to the 
presence of the emergent gauge field. The universal number that determines the DC conductivity is 
S/(0) = 0.010. It can be compared with Ej ^(0) = 0.085 in the absence of the gauge field, i.e. for 
the pure O(N) model. Extrapolating to the case of physical interest, N — 2, the conductivity reads: 

e 2 

cr b (0) = -x 0.020, 0(2) + damped gauge field (41) 



h 



x 0.170, pure 0(2) (42) 



We note en passant that this last number for the pure 0(2) model, 0.170, is very close to what 
was obtained in the small-e expansionPS 0.1650. Both these numbers lie near the self-dual value 
1/2tt ~ 0.159 which is associated with a conductivity equal to the quantum of conductance, e 2 /h. 

As the quantity that has a universal jump at the MIT is the resistivity and not the conductivity 
(unless we restrict ourselves to clean systems), we here give the expression for the rotor DC resistivity, 
Pb = l/<rb, 

p b = ~R, R = 49.8 (43) 
e 1 

or Pb — \ x 7.93 = 205 kO. This constitutes one of our main results: the value of the universal 
resistivity jump estimated in the large- N approximation, as shown in Fig. [TJ 

Although the exact numbers can only be trusted in the large- N limit, we expect that some semi- 
quantitative features are captured in our extension to N = 2. First, it is clear that the damped gauge 
field will necessarily make the conductivity smaller because it adds an additional scattering channel. 
The decrease in the rotor conductivity will lead to an increase in the resistivity, which translates 
into a bigger jump as one approaches the QCP from the metallic side, compared to a treatment 
that neglects gauge fluctuations. In our calculation, the increase is by a factor of ~ 8, which is 
substantial. Although, the actual enhancement might not be as large, our result suggests that the 
gauge fluctuations are the dominant source of current dissipation. Moreover, the strong scattering by 
the gauge bosons leads to the frequency dependence found in Fig. |8(b)| where the Drude-like peak 
occurring in the pure O(N) model disappears. Indeed, the small frequency conductivity is smaller 
than in the u/T>l limit. This "inverse Drude peak" might naively suggest that vortex excitations, 
whose conductivity is the quasiparticle resistivity, cr™ rtcx = p b = l/cr b: would be better suited, at 
least to describe electric transport. However, the vortices are known to give unreliable perturbative 
results for the supcrfluid-insulator transition in the pure 0(2) model. As the presence of the damped 
gauge boson does not alter the thermodynamic universality class, we suspect that this remains true 
in our model and that the dual vortex formulation does not offer any numerical advantage. 

We briefly comment on the precise frequency dependence of the conductivity at low frequency, as 
shown in Fig. [7] It is possible to determine almost exactly the analytic form of this frequency depen- 
dence for both the pure and gauged O(N) models, 29 which can be surprising given the complicated 
form of the QBE. Using the analytic expression, one can evaluate a low-frequency sum rule for the 
real part of the conductivity: J Q du a' b (uj) = constant. It is found that this integral equals the weight 
of the delta-function Drude peak obtained in the pure O(N) NLcM at N — oo, even in the presence of 
the gauge field. In particular this means that the inclusion of the interactions at 1/N (A or gauge bo- 
son mediated) only spreads the spectral weight of the delta function over a finite range of frequencies, 
this range being broader in the presence of the gauge field, as can be seen in the inset of Fig. [7] 



V. CONDUCTIVITY IN THE ENTIRE QC REGION 



We now examine how the conductivity changes as the system is tuned away from g — g c within 
the QC region, with g cx U/t, the ratio of the Hubbard repulsion to the electronic bandwidth. The 
QC regions is defined by T > A±, where A± are the two energy scales that vanish at the QCP. A + , 
defined for g > g c , is the Mott gap of the bosons, while A_, g < g c , is the phase stiffness of the rotors 
in their condensed phase. These two scales vanish approaching the QCP according to the power law 



(44) 
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FIG. 7. Universal scaling function for the real part of the conductivity above the QCP (S — 0). Q = Nuj/T 
is the rescaled frequency. The dotted line shows the scaling function in the absence of the gauge field. The 
inset shows for a larger range of w. 




(a) AT > 1 (b)N = 2 



FIG. 8. Frequency and temperature dependent rotor conductivity (real part), a) Sketch of a' b in the large- N 
limit. There is a Drude-like peak of width 1/N and height N. While it is due to incoherent scattering of 
thermally excited quasiparticles, the large-frequency 0(1) conductivity is from coherent and elastic scattering 
of field-excited carriers, b) Large- N result extrapolated to N = 2, the case of physical interest for the Mott 
transition. The solid blue (dotted purple) line corresponds to the conductivity with (without) the emergent 
gauge field. The dashed line at intermediate frequencies is a sketch of the expected crossover to the high 
frequency regime, where the gauge field becomes unimportant. 



where we are again using the signed energy scale associated with tuning the non-thermal parameter: 

5 = g- 1 -g- 1 ott/U-(t/U) c . (45) 

The dynamical exponent, z, is unity for all N, while the correlation length exponent, i>, depends on 
N: in the N — > oo limit, v — l/(d— 1) so that in d — 2, v = 1. In the QC region, the effective mass of 
the rotors, the saddle-point value of the A field, will change as 6 is varied. This will naturally affect 
the conductivity: as one approaches the SL, the effective mass of the charge excitations increases and 

this leads to a larger electric resistivity. The mass will depend on the ratio of A±/T: ^ = X± (t^) • 

At N — oo, we simply have A± cx |<5|, with the proportionality constants: 

A+ = 4tt<S, g > g c (46) 
A- = -5, g< g c (47) 

There, an analytic solution can be obtained for the mass scaling functiorPSl 

?-*(?)• < 48 » 

/ p 2wS \ 

X(6) = 2smh- 1 — , (49) 



which is plotted in Fig. [TU] 
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FIG. 9. Behaviour of the low temperature DC resistivity near the quantum critical (QC) Mott transition 
as obtained from the solution of the quantum Boltzmann equation. Panel c) shows the resistivity vs T for 
different ratios of the onsite repulsion over the bandwidth (tuned by 5). The corresponding cuts are shown in 
the phase diagram in a) and correspond to 5 — —0.01,-0.001,0,0.001,0.01 going from curve 1 to 5. Panel d) 
shows the resistivity vs S at different temperatures, with the corresponding cuts shown in the phase diagram 
in b). Curves 1,2,3 correspond to T = 0.5, 1.0, 2.5, respectively. The universal rotor resistivity at criticality is 
Pb = Rli/e 2 , with R = 49.8. (6 and T are given in a common and arbitrary unit of energy.) 



A. Results 



We have solved for the frequency dependent conductivity at different values of 6/T. The latter 
affects the polarization functions Id;, and n J b via the rotor effective mass. The polarization functions 
determine the propagators of the A and gauge boson, and thus the associated density of states the 
rotor excitations can scatter into. The universal scaling function for the low-frequency conductivity, 
£r, was numerically obtained in the QC region: 

e 2 fNui S\ 
a b =-NE^-,-y (50) 

In Fig. § we show the behaviour of the corresponding DC resistivity extrapolated to N = 2: pb(0) = 
l/<7b(0) = (fi/e 2 )/2£/(0, 5 zy /T). These numerical results should be compared with Fig. |3j where a 
sketch of the resistivity near the QCP was given, not restricted to the QC fan. It should also be 
compared with Fig. |15[ in Appendix |Aj showing the rotor contribution without the gauge field, as 



relevant for the conventional superfluid-insulator transition. Let us first examine Fig. 9(c) which 



shows the resistivity along constant-^ cuts (Fig. 9(a)). These curves can be naturally compared with 



sheet resistivity versus pressure experimental data, for instance, where S ~ P—P c , plays the role of the 
deviation from the quantum critical pressure. As mentioned above, the rotor resistivity corresponds 
to the total electronic resistivity relative to the residual value in the FL: pb = p — p m . Curve 3 shows 
the resistivity at the critical pressure (6 = 0): it is constant at low temperatures and takes a universal 
value AR, with R = 49.8. As one goes down in temperature at pressures differing from P c , the 
resistivity decreases approaching the FL (curves 4-5), or increases near the SL Mott insulator (curves 
1-2). The more pronounced resistivity jump p b — > occurs upon exiting the critical Fermi surface 
state and going to the (marginal) FL, as is shown in Fig. 3(c) Figure [9 (d) | presents the results from a 



complementary perspective: by fixing T and tuning pressure following cuts shown in Fig. |9(b)| This 
illustrates how the T = resistivity jump becomes smooth at finite temperature. By virtue of the 
scaling nature of the resistivity data, all curves cross at 6 = 0, where the resistivity is universal, ^R. 



16 




The resistivity data was extracted from a single universal scaling function 

Pb {u = 0) = —G Vj • ( 51 ) 

We emphasize that, in experiments, it is the resistivity relative to its residual metallic value should 
accordingly be examined for scaling in the vicinity of a quantum critical Mott transition. Notes: 
In obtaining the results shown in Fig. |9j we have replaced the N = oo scaling, S/T, by the one 
appropriate for N = 2: S zl/ /T, with v = 0.67 and z = 1 VW. (When 6 < 0, it is understood that 
5 ZV = sgn(<5)|<5| 2iy .) A caveat with the extrapolation is that we have used the mass scaling function 
obtained at N = oo. Although the specific form of m/T — X[8 ZV /T) will be different for N = 2, we 
expect it to be qualitatively similar, at least near criticality. 



VI. CONDUCTIVITY AT LARGE FREQUENCIES AND TEMPERATURES 

A. Large frequencies 

In this section we discuss the behaviour of the rotor conductivity in the large frequency limit 
uj m ~ T '. As discussed above, in the large- N limit, the finite universal conductivity in that region 
mainly results from the elastic, coherent transport of charged excitations created by the external 
field, as opposed to the incoherent transport of thermally excited quasiparticles relevant at small 
frequencies. In our above treatment of the quantum Boltzmann equation we cannot obtain that part 
of the conductivity as we have neglected precisely those processes that are dominant for w > T. 
Rather, the large-frequency conductivity a^jj — (e 2 /h)Y,ji(uj/T) can be obtained from a T = 
calculation, i.e. in the limit w/T — > oo. At N = oo, the A and gauge fields do not contribute and we 
recover the pure rotor contribution known from previous works! 2 -^ 



OJ 



1 



E' u (u>/T » 1) = ~ ^-*= °- 0625 (52) 
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where S d = 2/[r(d/2)(47r) d ' /2 ]. Contrary to the DC conductivity, which was infinite at N = oo, the 
fact that the high frequency conductivity is already finite in the free limit testifies about the different 
mechanisms at play, namely, its collisionless nature in contrast to the hydrodynamic transport at 
small frequencies. Including the 1/N correction that arises due to interactions of the rotors with the 
A field yieldsP 

Z' u (u>/T » 1) = 1 - |^ ^> 0.03998 (53) 

where r\ — 7V -1 (8/37r 2 ) is the leading correction to the anomalous dimension of the rotor field. We 
do not include the gauge fluctuations because they are made "massive" by the Landau damping. At 
T = 0, even the static component is ineffective, lying within the continuum of dynamical excitations. 
For a more precise estimate of E^ / (oo), we quote the Monte Carlo results of Cha et al.P^ Ej J (oo) = 
0.0454. This is the number we use for the large frequency conductivity in Fig. |8(b)"| 



B. Temperature dependence of electrical transport at criticality 



In our analysis of the quantum critical resistivity, we have so far limited our discussion to the 
universal quantum critical resistivity of the bosons. As the temperature is decreased to zero right 
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at the quantum critical point how is this universal value approached? To address this we first note 
that the full resistivity at low T is given by the Ioffe-Larkin formula as the sum of boson and fermion 
resistivities. At low-T the dependence of the fermion resistivity is dominated by spinons scattering 
off the gauge fluctuations. As mentioned in Ref. this leads to 




To obtain this dependence the gauge propagator at criticality, D 1 ~ + q, was used. Re- 

garding the subleading temperature dependence of the rotor conductivity, we have indications (see 
next subsection) that the treatment of the QBE requires more care, as in the case of a Fermi surface 
of spinons coupled to an emergent gauge fieldPDES We leave such analysis for future investigation. 
Notwithstanding, the leading non-constant T-dependence of the full resistivity will show a departure 
from the usual T 2 behavior because the rotor contribution is not expected to cancel the non-FL term 
provided by the spinons. This departure from a T 2 -dependence constitutes a further signature of the 
non-FL nature of the critical Fermi surface state. We note that the logarithmic correction due to the 
spinons will persist in the MFL phase, because the gauge field only becomes Higgsed by the rotor 
condensate at lower temperatures, where the usual FL is recovered. 



1. Necessity of careful treatment of subleading T -dependence of rotor conductivity 

To estimate the temperature dependence of the rotor conductivity at 5 = beyond its universal 
constant value, one can try to compute the transport scattering rate, just as was done for the spinon 



contribution, Eq. (54). We need to evaluate the imaginary part of the rotor self-energy due to the 
gauge fluctuations, with the usual additional factor of 1 — cos 9 in the integrand. The leading order 
dependence is universal and linear in T, in agreement with our QBE calculation. The subleading term 
is negative and goes like — T 2 //z, where /i is the fermionic chemical potential. We believe that such 
a negative contribution points to the inability of such a simple approach to capture the correct In- 
dependence. Indeed, this scattering rate would imply a conductivity that increases with temperature 
via <Tb ~ Tr t b r ~ 1/(1 — T/fi). 

Instead of using this semiclassical approach, we can turn to the full quantum Boltzmann equation. 
However, a straight-forward approach fails because the QBE contains terms that diverge due to 
scattering of low energy gapless gauge bosons. This can be related to the divergence of the rotor 
self-energy at finite T, which displays a logarithmic singularity in the infrared. It cannot be naively 
cutoff because the low energy properties of the gauge field and rotors were properly treated. This 
situation is exactly analogous to the case of fermionic spinons or non-relativistic bosons coupled to an 
emergent gauge field studied in the context of U(l) spin liquids.^ It was found that the divergence of 
the self-energy (or vanishing of the Green's function) is a natural consequence of the gauge fluctuations 
on a non-gauge-invariant operator, which in this case is the rotor Green's function. 

This type of singularity in the QBE was previously encountered in the context of the spinon-gauge- 
field problem by Kim et al. revisited latter by Nave and LeeP^Kim et al. found that a separation 
of the gauge fluctuations into static (with frequency less than the temperature) and dynamic ones 
together with the use of a gauge invariant momentum remove the singularities and allow the extraction 
of physical quantities. It was also noted that the static gauge fluctuations act as a random magnetic 
field that contributes to the conductivity by increasing the scattering of spinons, a situation very 
analogous to the effects of the gauge field on the critical rotors analyzed in this work. However, in 
the present system, the static gauge fluctuations dominate the low temperature transport and we 
did not need to take into account the dynamical or quantum gauge fluctuations, at least for the low 
temperature transport. The latter are essential to determine the higher temperature behaviour and 
one needs to perform a careful treatment analogous to Ref. |32j a task we leave for the future. 



VII. THERMAL CONDUCTIVITY 



The thermal conductivity also bears a signature of the critical Fermi surface, in which the emergent 
gauge field plays an even more important role than for electric transport. According to the Ioffe-Larkin 
composition rule, the thermal conductivities of the spinons and rotors addP^ k = Kb + k/. In other 
words, the thermal resistivities add in parallel, which is in contrast to the rule for electrical resistivity 
(addition in series) because of the absence of charge flow in response to the thermal gradient. In this 
sense, the thermal current is oblivious to the slave-rotor constraint relating the rotor charge to the 
spinon number, and the relative flow of the partons can proceed unconstrained. 

Let us first consider the rotor thermal conductivity, Kf,- In the absence of the gauge field, the rotors 
decouple from the spinons and their action reduces to the critical theory of the pure O(N) model. This 
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is a conformal field theory (CFT). On s ymm etry grounds, a CFT has infinite thermal conductivity 
because the energy current is conserved! 34 * 35 * Indeed, by virtue of conformal invariance, there exists 
a conserved and symmetric energy-momentum tensor: d^T^ — and T M „ = T UfJl , where /z, v are 
space-time indices. These two conditions imply that 4- J dx Ti T = 0, with spatial indices i = x,y. 
Put in words, the energy current, T ir , is conserved and will not be dissipated by interactions, contrary 
to the charge current. In our critical theory, this situation is avoided because of the damped gauge 
fluctuations, which naturally break conformal invariance and lead to a finite and universal critical 
thermal conductivity: 

k 2 

Kb = —jr-KT , (55) 

where if is a dimensionless number associated with the Mott QCP, just like R. The gauge fluctu- 
ations are more detrimental in the determination of the thermal conductivity than for the electric 
conductivity: whereas the later was already a finite universal constant without the gauge bosons, the 
former is formally infinite without the gauge fluctuations. In reality, this would not be the case due to 
the presence of irrelevant (in the RG sense) Umklapp scattering by the lattice, which would lead to a 
large, non-universal but finite conductivity. In the case under consideration, we do not need to refer 



to such processes because the gauge scattering is stronger and leads to a universal answer, Eq. (55). 

Although the full calculation of is beyond the scope of this work, we mention some of the 
important aspects. First, in the electric resistivity calculation performed above, a key simplification 
in the large-TV framework is that we can neglect the effect of the electric field on the A. As was 
explained in section |IV A| only one rotor flavor is directly coupled to the electric field so that its 
effects on the rotor polarization functions, which are obtained by summing over all rotor flavors, 
are subleading in N. This is no longer true when a thermal gradient is present, as all rotor flavors 
inexorably transport energy/entropy in the same way. Hence, the non-equilibrium corrections to the 
polarization functions cannot be neglected. To obtain the correct QBE describing the heat transport, 
one should use the Keldysh formalism, a task beyond the scope of the present work. 

We now turn to the spinon conductivity Kf. In the presence of weak disorder, the low temperature 
spinon thermal conductivity will scale like Kf = const x T, a form valid on both sides of the transition. 



Approaching from the FL, the constant k/T is simply La m , where L = ^- ( — ) is the usual Lorentz 



number while a m is the residual metallic conductivity, by virtue of the Wiedemann-Franz law obeyed 

k 2 

in the FL. As the critical point is reached, k/T jumps by a universal amount -§-K. Note that 
contrary to the electric conductivity, the thermal one is finite on the SL side, and is dominated by the 
spinon-gauge-field sector. 



A. Violation of Wiedemann-Franz 



Combining the electric resistivity and k/T jumps, we predict the Lorenz number, L = n/Ta = pn/T, 
will also jump at the transition, indicating a violation of the Wiedemann-Franz law. Approaching 

from the FL side, the Lorenz number will obey the Wiedemann-Franz law L = ^- (^f-) until it jumps 

by a universal amount (^f- ) 2 KR directly at the transition. 



VIII. DISCUSSION 



A. Experiments 



A preliminary analysis of unpublished pressure and te mp erature dependent resistivity data on both 
rt-(BEDT-TTF) 2 Cu 2 (CN);P and EtMe 3 Sb[Pd(dmit) 2 ]j221 provide encouraging hints regarding the 
presence of a quantum critical Mott transition as described in this work. Indeed, at sufficiently high 
temperatures, the sheet resistivity saturates to a constant ~ h/e 2 , as we predict. The qualitative 
pressure and temperature dependence is also roughly in agreement with our results. A closer analysis 
of the data will be needed to make a stronger statement. Not to mention the potential complications 
with the inhomogeneous effects of pressure on the relatively soft organic salts. 



B. Disorder 



We discuss the effects of disorder and provide an estimate for the temperature at which we expect 
it to become important in the context of the organics. We shall assume the disorder is weak as is 
appropriate for these materials. Then its main effect will be to modify the gauge field's inverse Green's 
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function to (we work in imaginary time) 

o y^T^v (56) 

y/q 2 + l/P 

Here, I is the spinon elastic mean free path, and we have reinstated c, the rotor velocity, which is 
on the order of Vp. er° is a shorthand for the T = rotor conductivity, <7b(w/T 1). The zero 
Matsubara frequency gauge mode does not respond to the change in the Landau damping, which 
vanishes at zero frequency. So the question is about the effect on the non-zero Matsubara frequency 
components. These will notice the finite mean free path when q < l/l or relating the gauge boson 

energy (~ T) to the momentum via the z = 2 scaling Q ~ q 2 , we get T < jprj. Using c f» Vp, this 
can be cast as 

<■"» 

where /i is the Fermi energy or chemical potential. The organics are good metals, so we take kpl ~ 10, 
although the actual mean free path is probably larger. Also, [i w 10 3 K and Ub ~ 0.05, in units of e 2 /h. 
We thus conclude that within our framework, using parameters relevant to the organics of interest, 
disorder will start playing an important role at temperatures 

T<0.5K. (58) 

Above that temperature, we can treat the higher Matsubara modes without including disorder, and 
the results of this paper should provide a valid description. 



IX. CONCLUSION 



We have analyzed the transport signatures of a quantum critical Mott transition between a Fermi 
liquid metal and a paramagnetic Mott insulator in two spatial dimensions, at fixed filling. In this 
scenario, the Mott phase is characterized by a Fermi surface of spin-only quasiparticles, resulting in 
a gapless U(l) spin liquid. The physics of such a transition is conveniently captured by a slave- rotor 
field theory of the electronic Hubbard model, in which the charge fluctuations are described by a 
quantum XY model of rotors, coupled to an emergent gauge field. The superfluid phase of the rotors 
corresponds to the metal, while a Mott insulator results in the disordered phase, where only spin 
fluctuations remain gapless. Directly at the transition, a strongly correlated non- Fermi liquid metal 
emerges where the electronic Fermi surface is on the brink of disappearance; it is an instance of a 
critical Fermi surface. The zero temperature electric resistivity of such a critical state was predicted 
to be greater than that of the Fermi liquid by a universal amount, Rh/e 2 , where R is a dimensionless 
number associated with the Mott quantum critical point, see Fig.[T] We have obtained an estimate for 
this number, R « 49.8, via the solution of a quantum Boltzmann equation for the charge fluctuations, 
analyzed in a large- TV limit. We found that the emergent gauge fluctuations strongly contribute to 
the universal resistivity jump, albeit being ineffective at changing the universality class of the charge 
fluctuations from 3D XY. This is so because the static gauge fluctuations escape the Landau damping, 
and are responsible for the strong, and universal, scattering of the critical charge fluctuations. Though 
we have focused on transport properties which are best accessed in experiments, we may anticipate 
that the static gauge field fluctuations might also affect other dynamical quantities (for instance the 
low frequency form of the electron self energy at non-zero temperature) in the quantum critical regime. 

We have further examined how this resistivity jump evolves at finite temperature and as one 
changes the ratio of the Hubbard repulsion to the electronic bandwidth (experimentally tunable 
by applying pressure). We have obtained a universal scaling function that can be used to col- 
lapse the temperature and pressure dependent resistivity in the quantum critical Fermi surface state: 
p — Pm — {h/e 2 )G{8 zu /T), where p m is the residual resistivity in the Fermi liquid, 6 can be mapped to 
the deviation from the critical pressure, and z and v are the usual dynamical and correlation length 
exponents, respectively, of the 3D XY universality class. In particular, G(0) = R. 

Turning to thermal transport, we make the prediction that the low temperature thermal conduc- 
tivity (divided by temperature), k/T, has a universal jump at the transition, Fig. [I] The gauge 
fluctuations are responsible for this jump by breaking the conformal invariance of the charge fluctu- 
ations. Together with the jump of the electric resistivity, the thermal conductivity jump leads to a 
violation of the Wiedemann-Franz law by the critical Fermi surface state. 

Regarding experiments, the organic salts k-(BEDT-TTF) 2 Cu 2 (CN)3 and EtMe 3 Sb[Pd(dmit) 2 ]2 are 
candidate materials for the transition covered in this work. Indeed, at ambient pressure they display 
numerous signatures characteristic of a gapless quantum spin liquid, and become metallic under the 
application of hydrostatic pressure. As mentioned in section |VIII| recent unpublished data for the 
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temperature and pressure dependent resistivity near the Mott transition seem to indicate a jump on 
the order of h/e 2 near the transition, as well as qualitative agreement with the scaling we propose. A 
confirmation would point towards the first experimental example of a quantum critical Mott transition 
(of fermions), as well as a non-trivial signature of spin-charge separation in two dimensions. 

On the theoretical side, it would be interesting to investigate the charge transport within a controlled 
dimensional expansion, with small parameter 3 — d where d is the spatial dimension, and see how the 
results compare with the large- TV expansion. Numerical simulations directly in d — 2 and at N = 2 
would also be desired. Regarding the jump of the thermal conductivity, a full treatment within the 
large- TV expansion is still missing, due to the complication of drag of the constraint field of the 0(N) 
non-linear sigma model. A first step would be to establish the divergence of the thermal conductivity 
in the pure XY model, as required by conformal invariance, then to include the emergent gauge 
fluctuations. Taking a broader perspective, we envision that the results regarding the effects of the 
damped gauge fluctuations on the transport of relativistic or particle-hole symmetric quasiparticles 
discussed in this work can be applied to other systems. 
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Appendix A: Transport in the pure O(N) rotor theory 

In this appendix we present details regarding transport properties of the pure 0(N) NLctM, where 
the N = 2 theory describes the superfluid-insulator transition of bosons in two dimensions, as present 
in the Bose-Hubbard model at integer filling for instance. Some of these results correct previous ones, 
or cannot be found in the literature. In addition, they are important as a comparison ground with 
the rotor theory used to describe charge criticality in the Mott transition under consideration. The 
latter is a gauged version of the O(N) model, with a Landau damped gauge boson. 

The bare action for the pure O(N) theory is: 

S b = ~J{\d v b a \ 2 + i\(\b a \ 2 -l)) , (Al) 

where are a runs from 1 to N/2, so that the N/2 complex fields form a real O(N) vector field. 
The above action has full O(N) global invariance. (The notational choice of using N/2 complex fields 
instead of N real ones is more convenient from the point of view of the slave-rotor formulation.) Again, 
A is a Lagrange multiplier field enforcing the constraint that the rotor is unimodular: ^ Q \b a \ = 1. In 
the N = oo limit, the field theory becomes free, allowing a 1/N perturbative expansion. At N = oo, 
the expectation value of the A field plays the role of an effective mass for the rotors in their insulating 



phase, as described in section III A The saddle point equation for A is still given by Eq. (12), because 
the fluctuating gauge field does not affect the N = oo rotors. In the symmetry broken phase, the 
saddle point equation for the effective mass has no solution, indicating long range order (at T = 
only for N > 2). 

We now turn to the charge transport near the quantum critical point. Out of the many charges 
present at large TV, we couple the electric field to a single one, say the complex component b\. We 
focus on the regime where the external field has a small driving frequency compared with temperature; 
this corresponds to the collision dominated hydrodynamic regime. At N = oo, the DC conductivity is 
infinite because interactions are entirely suppressed in that limit. At order 1/iV, the A field acquires 
a non-zero propagation amplitude and can scatter the positive and negative charge excitations, the 
holons and doublons, respectively. The quantum Boltzmann equation for the holons and doublons 



can be read off Eq. (24), where one simply has to set the gauge propagator, oc D(fl,q), to zero. As 



explained in the main text, we expand the holon and doublon distribution functions to linear order in 



21 



the electric field: f±(k,u>) = 2Tr8(oj)n(£ k ) ± E ■ kip(k,Lu). Linearizing the QBE, we obtain 

-iu<p{p,w) +g( P )/T 2 = ~ l-F x (p)<p(p,w) + f dp'K x {p,p')<p(p',Lj)\ (A2) 



A( P ,p') = 




dp' 


l\P + P'\) 


A c ( P ,p') = 


C dd An 




|,|p + p'|)(-cos0) 


Adh{p,p') = 


f 2 " dd A, 
J* 2^ 


+ £ P >, 


\P + P'\) 


A,dh(P,P') = 


f 2 " M A, 

Jo 27r 


+ £p', 


\p + p'\)(-cos9) 



which is the same as Eq. (31), except that we have set the scattering rate due to gauge fluctuations, 
F a , to zero. As before, g(p) = d ep n(e p )/e p . The functions F x and K x can be written in the simple 
form: 

r°° v'dv' 1 

Fx(p)= V -^- [AfartfrM^+AtoipJhtoM)] (A3) 

ill CpCp' 

KxM) = f— [AMhip 1 ',P)\- A cA M) lAh { P > »] t (A4) 
where we have defined the two 7 functions: 

1 - e~V 

7(p,p ) = n(e p , - e p ) - n{e p .) = ^ _ _ (A5) 

7dh(p,p') = n(e p i) - n(e p + e p >) (A6) 
and the four A functions: 

(A7) 
(A8) 
(A9) 
(A10) 

where 9 is the angle between p and p' . We have also defined the scaled spectral function of the A 
field: A(Q,q) = — 3?[l/rLj(Q, q)], where the factor of 2/N in the A propagator, 2/NHb, was omitted 
from the definition of A. The sign of A is such that A(ft > 0, q) > 0. The A functions can be readily 
seen to be positive and symmetric. The "c" indicates a cosine in the angle integral, while "dh" stands 
for doublon-holon. To understand this last piece of notation, we need to examine lib, specifically its 
imaginary part, 

n* («.«) = / i 6 jf fc g fc+g (K £ fc+q) - n(e k )\S(\e k+g - e k \ - O) 

+ (1 + n(ejfe) + n(e fe+9 ))5(e fc + e k+q -n)}, (Aii) 

which is non-zero in two separate regions, a low and high energy one: Q < q and 57 > ^/ (2m) 2 + q 2 , 
respectively, as is illustrated in Fig. [TT] The latter region arises from the production of holon-doublon 
pairs, which requires an energy beyond a certain threshold. The A-boson spectral function has support 
in the same regions, being proportional to n£: A = -$[1/II&] = U^/[(U' b ) 2 + (U' b r ) 2 }. We plot the 
numerically-evaluated spectral function in Fig. [l2j We note that A vanishes at the boundaries of the 
region where II' b ' = 0. Going back to the A functions, it can be seen that those labelled "dh" receive 
contributions precisely from the doublon-holon pair production region. 
By performing the rescaling Co = toN/T and 

S(p,w) = ^>(P,w) (A12) 
we obtain a universal, parameter-free equation: 

-iw*(p,w) +g(p) = -F x (p)<S>{p 7 £o) + J dp' ' K x {p,p')${p' ' ,u) (A13) 

We solve the integral equation numerically by discretizing the momentum variables and expanding 
the unknown function $ in terms of Chebyshev polynomials. This procedure converts the integral 
equation into a matrix equation, which is solved by simply inverting the matrix corresponding to 



the kernel K\. We plot the solution in Fig. 13 3>(p, ui) goes to a constant as p — » and decays 
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FIG. 12. a) Spectral function for the A-bosons: A(Q,q) — — O^l/II^fi, q)]. b) The corresponding density 
plot. 



exponentially for p> 1, i.e. for momenta much greater than the temperature. Also, as the driving 
frequency ukw goes to zero, the imaginary part, vanishes as expected since the DC conductivity 
is purely real. 

The resulting frequency-dependent conductivity can be obtained by integrating <I>(p, Co) over mo- 
mentum, as shown in Eq. (39 1. As a result we obtain a scaling function for the co nduc tivity: 
<7b = (e 2 /h)N^i(ui, 8 ZU /T). The numerical solution for the scaling function is shown in Fig. [l4j where 
a) shows the frequency dependence of both the real and imaginary parts of the scaled conductivity 
for <5 = 0, while b) shows the real part away from criticality at five different values of 5/T. Fig. |14(b)| 
shows that the conductivity increases going away from the disordered to the symmetry-broken phase, 
corresponding to sweeping 8 from negative to positive values. 

For N = 2, the conductivity is given by ct& = (e 2 /fi)2E/(2o;/T, S zu /T), so that the curves in Fig 
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give the conductivity (in units of 2e 2 /h) versus 2 x uj/T. In general, the scaled frequency Co on the 
horizontal axis should be multiplied by 1/N to get the dependence on lo/T. In the large- A limit, we 
recover a Drude-like peak, as the charged quasiparticles interact more weakly. As mentionned at the 
end of Section |IV B the essentially exact frequency dependence of the small frequency conductivity, 



cf. Fig. 14(a) was recently determined!^ This result can be used to prove a sum rule constraining the 



integral of the real part of the conductivity to the weight of the delta- function Drude peak at N = oo, 
i.e. in the free limit of the theory. This provides an excellent check on the numerics. 



Focusing on N — 2, Fig. 15 shows the behaviour of the resistivity in the QC region. We have 



indicated the presence of a Kosterlitz-Thouless transition by a solid line at the interface between the 
QC regime and the SF. For N > 2, this finite- T phase transition is converted to a crossover. This 
figure should be compared with Fig.[9j where there the bosons are coupled to a Landau damped gauge 
field. We see that not only does the gauge field make the resistivity larger, but its variations are also 
more pronounced as S and T are changed. 



Appendix B: Thermal conductivity 



We discuss some aspects of the thermal conductivity of the rotors. In particular, via the solution of 
a QBE, we explicitly show that in the absence of gauge fluctuations, the rotors have infinite thermal 
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FIG. 13. Solution for electric-field induced deviation to the rotor distribution function, $(p, uj), for ti) = 
Nuj/T — 0.5. The solid (dashed) line shows the real (imaginary) part, multiplied by p 2 . 




(a) (b) 

FIG. 14. Conductivity universal scaling function for the pure O(N) theory, a) Frequency-dependent scaling 
at finite T above the QCP, 8 = 0: Ob = NF<i(uj), where the argument is uj — Nuj/T. The solid (dashed) line 
corresponds to the real (imaginary) part, b) The scaling function (real part) as a function of the departure 
from criticality, S/T. 



conductivity, as expected based on conformal invariance. We then go on to argue that with the 
addition of the damped gauge field, the conductivity becomes finite, specifically «;& = (k%/H)KT, 
where if is a dimensionless constant. In the large- N formulation one needs to include effects of the 
thermal gradient on the constraint field A, complicating the analysis. We thus consider an alternate 
but equivalent formulation of the critical theory where the hard constraint of the XY rotor model, 
\b\ = 1, is traded for a soft one via the addition of quartic term |6| 4 . 



1. Thermal response of pure 0(2) model: CFT and zero modes 



The critical theory for the "soft" rotor model, without gauge field, is described by the action 

S= 1 -Jd^ X \d^+m^ + ^\b\\ (Bl) 

where d = 3 — e is the space-time dimension, and to 2 = e(47r 2 /15)T 2 and u = e(247r 2 /5). This is the 
standard Wilson-Fisher fixed point resulting from a perturbative RG treatment (see for e.g. [2D which 
studies the electric transport properties of that model). We have thus traded the large- N expansion 
for a dimensional expansion in e, N = 2 being fixed. 

We are interested in the linear response of the bosons to a thermal gradient. As we did above in the 
case of the electric transport, we consider the QBE for the holon and doublon distribution functions, 
f±. The static thermal gradient is imposed by including a position dependent temperature, T(x), 
which in turn generates position dependence for the distribution functions f±. The corresponding 
QBE, which can be found in RefsI3"51 and EM reads: 

« & ~=I± [/+,/-], (B2) 



where Vk = dek/dk and = \Jk 2 + to 2 . The collision term is different from the large- TV expansion 
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FIG. 15. Behaviour of the low temperature DC resistivity near the quantum critical (QC) superfiuid- 
Mott-insulator phase transition. In a) and b), the solid line delimiting the SF indicates a Kosterlitz-Thouless 
transition. Panel c) shows the resistivity vs T for different ratios of the non-thermal parameter 8, which would 
tune the onsite repulsion, U , over the tunneling amplitude, t, in a Bose-Hubbard model, for example. The 
corresponding cuts are shown in the phase diagram in a) and correspond to S = —0.01,-0.001,0,0.001,0.01 
going from curve 1 to 5. Panel d) shows the resistivity vs S at different temperatures, with the corresponding 
cuts shown in the phase diagram in b). Curves 1,2,3 correspond to T = 0.5, 1.0, 2.5, respectively. The universal 
resistivity at criticality is pb = Rh/e 2 , with R = 5.88. (5 and T are given in a common and arbitrary unit of 
energy.) 



as it arises from the quartic interaction: 

o„2 rJL -72 ju. C2ir) d+1 
I± = — g- / II (^(^ - ^±) L -^S(k + fei - fe 2 - k 3 )S(e + e kl e k2 - e ka ) (B3) 

where scattering processes out of state {fc, =b} have 

?T = 2/±(fc)/ T (fci)[l + /±(fca)][l + Mk 3 )] + /±(fe)/±(fei)[l + /±(fc 2 )][l + f±(k 3 )} , (B4) 

J 7 ™ is obtained by interchanging f± and 1 + f±. We now linearize the QBE to linear order in the 
temperature gradient. From the LHS, 

-g^ = -e k de k n(e k )— , (B5) 
while the expanded distribution function reads 

/±(fc) = n(e fc ) + (B6) 

with <fi characterizing the departure from equilibrium due to the applied thermal gradient; it is analo- 
gous to ip used above in the context of electric transport. An important difference is that the thermal 
gradient leads to the same non-equilibrium distribution functions for both the positive and negative 
charge excitations: /+ = /_, contrary to the case with an electric field. In the linear-response regime, 
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the heat current is simply the energy current: 

Jh = I v k e k [f + (k) + f_(k)] = 2~ / ; (B7) 

Jk 1 Jk 

while we naturally get a vanishing electric current because the latter involves the integral of f+(k) — 
f-(k), see Eq. (35 1. This is as it should be since the thermal conductivity is measured in an open 
circuit setup in which there is no electric charge flow in the steady state. 

The linearized QBE reads 

-d k n k = I[ct>{k)] , (B8) 
where the functional for the linearized collision integral is 

I[4>(k)] = J dk 1 dk 2 Ii(k,ki,k 2 )n kl n k2 (l + n k2+ k 1 -k) (B9) 

-6(1 + n k ) / dk 1 dk 2 — -h{k, fci, k 2 )n k2 n k+kl - k2 (BIO) 
J n kl 

-12n k [dkidk 2 — l3(k,k 1 ,k 2 )n k2 (l + n k2+k - kl )X. (Bll) 
J n kl J 

The functions Ii result from the angle and £3 angle integrations of the (5-functions and are given in 
Ref. [2"T1 We are using the shorthand n k = n(k), having dropped the mass, 0(y/e), in the dispersion 
relation, e k = \Jk 2 + m 2 ps k, to leading order in e. 

By conformal invariance we expect the thermal conductivity, which is given by an integral over <p, 
to be infinite as explained in Section VII If we could invert the linear functional I [(/)], we would obtain 
4>(k) = d k n k ], yielding a finite function </>(fc), as is the case for the electric conductivity. This 

would in turn imply a finite thermal conductivity, running against general symmetry arguments. In 
fact, I[4>] cannot be inverted. In other words, the linearized scattering integral has a zero mode, i.e. a 
zero eigenvalue eigenfunction <^o(fc) such that /[</> (fc)] = 0. This zero mode turns out to be the LHS 
of the linearized QBE: 

<M fe ) = -d k n k = n k (l + n k ) , (B12) 

as we have verified by direct substitution. Hence, the equation has no solution, or formally <f> = 00, 
which implies that the CFT describing the Wilson-Fisher fixed point has infinite thermal conductivity. 
It is an excellent check on the formalism, as it is a priori not obvious how this divergence would come 
about within the QBE framework. 



2. Thermal response of gauged 0(2) model 



Let us now return to the case of relevance for the QC Mott transition. In that case the rotors are 
coupled to a damped gauge field. As was explained in the main body, only the static gauge fluctuations 
are effective at scattering the rotor excitations. These provide a simple elastic scattering term to the 
QBE: 

-d k n k = I[<t>(k)], I[<i>}=I[<p\-^- (B13) 

t q 

Although / itself is not invertible, I is for any finite value of the scattering time r a . In particular, the 
zero mode of I no longer is one for /: I[—d k n k ] = d k n k /T a , i.e. it now has eigenvalue — l/r . 



We have solved the QBE, Eq. (B13l, numerically and have verified the presence of a singularity 
when the gauge scattering rate is absent. (For simplicity, we have used a momentum- independent 
scattering rate but the same conclusions will hold generically.) More precisely, the numerics yield a 
very large value of k/T ~ 7 x 10 4 in the CFT limit, i.e. when l/r a = 0. This value grows as we 
increase the number of Chebyshev polynomials used to expand <f> suggesting a divergence in the limit 
where an infinite number of basis polynomials is used. Moreover, it becomes 0(1) even for very small 
values of l/r a , consistent with the fact that the gauge scattering rate moves the system away from 
the conformal point thus rendering the thermal conductivity finite. 
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3. Similarity with fermionic CFT of Dirac fermions 



A similar analysis of conformal zero modes in the context of thermal transport in a fermionic 
CFT of Dirac fermions was performed in Ref. 40, A zero mode was identified in the corresponding 
linearized QBE, the analogue of Eq. ( B8 ) . It was was found that the zero mode is essentially given 
by <fio{k) = —dkriF(k) = nf(fc) (1 — nF[k)), where np is the Fermi-Dirac distribution. This is the 
fermionic analogue of the bosonic mode obtained above, Eq. (B12|. It was further found that the 
introduction of anisotropy for the Dirac fermions breaks conformal invariance and makes the zero 
mode massive. As such, the anisotropy can be seen as the analogue of the gauge scattering rate in 
our case. 



Appendix C: Rotor current polarization function 



The static, non-regularized rotor current polarization function is obtained from: 

d 2 p (2q x p) 2 1 



{2nf 



6^ LxP 1 



(Cl) 



P+q 



where e p = \Jp 2 + m 2 is the rotor energy. Note the presence of a factor of 4 due to the (squared) 
vertex between the transverse component of the gauge field and the rotor current. Performing the 
Matsubara sum we obtain 



4(0,9) 



d 2 p 1 (2q x pf 



(2tt) 2 2 e 2 



p+q 



1 + 2n(e p+q ) _ 1 + 2n(e p ) 

e p+q e p 



(C2) 



This integral is UV-divergent; putting a cutoff would lead to n J 6 (0,0) 7^ 0, which would violate the 
U(l) gauge invariance. We regulate by subtracting the integrand evaluated at q = 0. Slight care must 
be used in doing so because of an undetermined limit 0/0. The regulated expression reads 



4(0,9) = 2 



d 2 p 

W) 2 



(q x p) 



-p+q 



1 + 2n{e p+q ) 1 + 2n{e p ) 



c p+q 



1 

2~c 



»[1 



2n(e p ) - 2e p n'(e p )] j 
(C3) 



where n'(e) — d e n(e). This integral can be evaluated numerically and the result is plotted in Fig. 16 
As expected, for q 3> T the polarization function scales like q, as in that regime the mass (~ T) 
is negligible and we recover the zero temperature behaviour. At small q, we obtain a q 2 scaling, 
consistent with the diamagnetic response of massive bosons. 

We can num eric ally extract the slope of the polarization function at large q/T to obtain the result 
quoted in Eq. (181, of ~ 0.063. This corresponds to the rotor conductivity in the uj/T ^> 1 limit, and 
agrees with the analytic expression 1/16 = 0.0625 given in Section VIA 
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